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Recent observations by the Swift satellite have revealed long-lasting 10 2 — 10 5 s), “plateau-like” X-ray 
afterglows in the vast majority of short gamma-ray bursts events. This has put forward the idea of a long-lived 
millisecond magnetar central engine being generated in a binary neutron star (BNS) merger and being respon¬ 
sible for the sustained energy injection over these timescales (“magnetar model”). We elaborate here on recent 
simulations that investigate the early evolution of such a merger remnant in general-relativistic magnetohydro¬ 
dynamics. These simulations reveal very different conditions than those usually assumed for dipole spin-down 
emission in the magnetar model. In particular, the surrounding of the newly formed NS is polluted by baryons 
due to a dense, highly magnetized and isotropic wind from the stellar surface that is induced by magnetic field 
amplification in the interior of the star. The timescales and luminosities of this wind are compatible with early 
X-ray afterglows, such as the “extended emission”. These isotropic winds are a generic feature of BNS merger 
remnants and thus represent an attractive alternative to current models of early X-ray afterglows. Further impli¬ 
cations to BNS mergers and short gamma-ray bursts are discussed. 


I. INTRODUCTION 

Coalescing binary neutron stars (BNSs) and neutron star- 
black hole (NS-BH) binaries are not only among the most 
promising candidates for direct detection of gravitational 
waves (GWs) with upcoming ground-based interferometers 
such as advanced LIGO and Virgo HI. The accretion of a 
remnant torus around a BH formed promptly after merger also 
represents the leading scenario to explain the formation of a 
short gamma-ray burst (SGRB; e.g., El). However, the re¬ 
cent discovery of long-lasting X-ray afterglows (referred to as 
“extended emission” and “X-ray plateaus”; e.g., El) in the 
vast majority of SGRB events observed by the Swift satel¬ 
lite |4) indicate ongoing energy injection on timescales of up 
to ~ 10 4 s, which are incompatible with the short accretion 
timescale of a remnant torus (< Is). 

These X-ray afterglows can, however, be interpreted by as¬ 
suming that the merger leads to the formation of a long-lived, 
rapidly rotating NS (a millisecond magnetar; henceforth re¬ 
ferred to as the “magnetar model”; e.g., 0 ). Despite indi¬ 
vidual differences in the phenomenology, different variants of 
the magnetar model have in common that they postulate the 
formation of a uniformly rotating, highly magnetized NS and 
that at least the late-time X-ray emission is powered by dipole 
spin-down radiation, which requires an essentially baryon- 
free environment around the NS. 

The magnetar model challenges the NS-BH progenitor 
channel for SGRBs, since a magnetar cannot be formed in this 
case, and we are thus led to consider BNS mergers as the main 
formation channel. The outcome of a BNS merger depends on 
the equation of state of nuclear matter at high densities, which 
is unknown. However, the recent discovery of high-mass neu¬ 
tron stars ID indicate a maximum mass for stable NSs of 
Mrov > 2 M 0 . Furthermore, the mass distribution of BNSs 
is peaked around 1.3 — 1.4 M 0 and this yields an approxi¬ 
mate mass of the merger remnant of ~ 2.4 M 0 Q. Combin¬ 
ing this with the fact that uniform rotation can centrifugally 
support NSs against gravitational collapse up to masses of 
~ 1.2 Mtov ^ 2.4 M 0 (5), merger remnants are very likely 


supramassive. If the merger remnant is only slightly hyper- 
massive (i.e. above the maximum mass for uniformly rotating 
NSs), it can still migrate to a supramassive star by mass loss 
through early winds while still supported by differential rota¬ 
tion (see below). Supramassive NSs are typically long-lived 
and the exact lifetime depends on the timescale to remove a 
significant fraction of the rotational energy (e.g., via dipole 
spin down). Hence, it is conceivable that in the vast majority 
of BNS mergers a long-lived NS is formed, which supports 
the magnetar model. 

From simulations it is known, however, that at early times 
the BNS merger remnant is characterized by a high degree of 
differential rotation and thus substantial magnetic field rear¬ 
rangement and amplification. These conditions are very dif¬ 
ferent from pure dipole spin down of a uniformly rotating 
magnetar, as assumed in most magnetar models. This leads 
us to investigate this early phase of a differentially rotating 
BNS merger remnant by fully general-relativistic magneto¬ 
hydrodynamic simulations. We elaborate here on recent re¬ 
sults of such simulations 0, which revealed the generation 
of a very dense, highly magnetized wind that carries sufficient 
Poynting flux to explain typical luminosities of early X-ray 
afterglows. The present paper focuses on discussing the wind 
properties in more detail and its implications for SGRBs and 
their X-ray afterglows. 


II. INITIAL DATA AND NUMERICAL SETUP 

Employing the general-relativistic ideal magnetohydrody¬ 
namics code WhiskyMHD fTOl combined with the publicly 
available Einstein Toolkit and its spacetime evolution code 
McLachlan ED, models of hypermassive neutron stars 
constructed with the RNS code E3 were evolved, which we 
endowed with different initial magnetic field geometries (for 
details, see (5)). These models are representative of the typical 
outcome of a BNS merger, once the merger remnant has set¬ 
tled down to a roughly axisymmetric state and assuming that 
the merger has not led to a prompt collapse to a BH. Here, we 
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FIG. 1. Magnetic field strength (left) and rest-mass density (right) of the wind from a differentially rotating NS (corresponding to model rand 
of [9|; the star is indicated by black density contours in the lower-left corners). These snapshots showing a very dense, highly isotropic and 
highly magnetized wind are taken at the end of the simulation. 


elaborate, in particular, on model rand of ED, which has a 
“random” initial magnetic field geometry inside and outside 
the star and which we consider to most closely resemble the 
actual magnetic field geometry resulting from a BNS merger. 
It is important to stress that our conclusions in this paper and 
that of (9) mostly depend on differential rotation and thus also 
apply to differentially rotating supramassive NSs (which to¬ 
gether with hypermassive NS are the most likely result of BNS 
mergers, see Section [I]). 


III. EVOLUTION, MASS EJECTION, 
ELECTROMAGNETIC EMISSION 

The early evolution of a differentially rotating BNS merger 
remnant is characterized by magnetic winding: poloidal mag¬ 
netic fields in the interior of the star are “wound up”, produc¬ 
ing a strong toroidal component and enhancing the magnetic 
pressure. Rotational energy is thus converted into magnetic 
energy. We find that within a few average rotational periods, 
the increase of magnetic pressure in the surface layers of the 
star is sufficient to overcome the gravitational binding of ma¬ 
terial in the vicinity of the stellar surface. This results in a 
highly magnetized, highly isotropic and very dense wind of 
matter from the NS that heavily pollutes the environment with 
baryons (cf. Figure [l} f] This environment of the NS is very 
different from an essentially baryon-fee pulsar wind that is 
usually hypothesized in magnetar models. 

A global measure of the baryon pollution is depicted in the 
rightmost panel of Figure[2j which shows the total mass cross¬ 
ing a coordinate sphere of radius r = 120 km and r = 150 km 


from the centre of the star per unit time during the evolution 
of the simulation. The wind at these distances is still rather 
turbulent during the later part of the evolution, which results 
in a noisy mass-loss rate. Nevertheless, the mass-loss rates 
as measured at different distances from the star roughly agree 
and yield a stationary mass-loss rate of a few in 10 _3 solar 
masses per second (M rs j 10 -3 M 0 s -1 ). 

The morphology of the wind is depicted in the left and mid¬ 
dle panel of Figure[2]in terms of the rest-mass density profiles. 
The directionally averaged profile^] at different times (e.g., at 
50 ms and 60 ms) agree very well up to the outer radius of the 
bulk of the ejected matter at that particular time. A fit to these 
profiles reveals a characteristic average wind density profile 
p oc r~ 2 that is constant in time. According to mass conser¬ 
vation, we have at fixed radius 


M = 


9M(r, t) 
dt 


Airr 2 p(r, t)v(r , t ), 


( 1 ) 


where M is the amount of mass inside the volume of ra¬ 
dius r and v(r,t ) is the radial velocity of the wind at time 
t and radius r. Combining p(r, t) = p(r ) oc r ~ 2 and 
M(r, t) = M = const thus yields u(r, t) = Vo = const, 
i.e. an overall bulk speed of the wind that is constant in time 
and space. Assuming M^5x 10 -3 M 0 (see above and Fig¬ 
ure [2]), we obtain a subrelativistic wind expansion speed for 
the model considered, Vo ~ 0.03 c, where c is the speed of 
light. 

Although the mass-loss rate is roughly constant during the 
timescale of up to 60 ms that we can access through our sim¬ 
ulations, it is expected to change significantly at later times. 
The timescale for a change of M and thus for a change of 


1 In the case of a dominant dipole moment in the initial magnetic field struc¬ 
ture, we expect an additional collimated wind component directed along 2 The system is very close to axisymmetry and averaging over the polar angle 
the rotation axis (9). is sufficient to obtain a directionally averaged profile. 
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FIG. 2. Left and middle panel: rest-mass density profiles of the wind from a differentially rotating NS (corresponding to model rand of 0) at 
different times towards the end of the evolution (density profiles along the x and z axes in cyan and purple, respectively, polar angle averaged 
profiles in green and blue). The black line shows a fit p oc r 2 to the angular averaged profiles. Right panel: instantaneous mass-loss rate 
associated to the wind as measured at different radii during the evolution. 


the wind properties depends on the timescale for removal of 
differential rotation (see Section [TV]). Furthermore, neutrinos 
radiated from the hot interior of the NS can heat and drive 
an outflow from the surface of the NS and thus lead to addi¬ 
tional mass loss (e.g., lH3l ). Depending on the magnetic field 
strength and the rotation period, however, the magnetically- 
driven wind discussed here may well be the dominant process 
for mass loss. 

The Poynting flux associated to the magnetically-driven 
wind in our simulations corresponds to stationary electromag¬ 
netic (EM) luminosities 

i EM =* 10 48 B'l- Rl & Pzl erg s" 1 , (2) 

where B is the average magnetic field strength in the outer 
layers of the star once the system has reached its stationary 
state, R e denotes the equatorial radius, and P is the (central) 
spin period m. This relation is universal in the sense that it 
does not depend on the initial magnetic field geometry and it 
is different from dipole spin-down emission considered in the 
magnetar model, in which case Lem oc B^ R® P -4 , where 
P p is the dipolar field strength at the pole. Assuming that 
only a small fraction rj = L^/Lem ~ 0.01 — 0.1 of the 
luminosity (Equation® is converted into X-rays, the observed 
X-ray luminosities L^| ^ 10 46 — 10 51 ergs -1 of extended 
emission and X-ray plateaus of SGRBs (e.g., C3i) are obtained 
by magnetic field strengths B 10 14 — 10 17 G, which are 
likely produced during and after a BNS merger [9 , 14 ]. 


IV. CONCLUSIONS: CONSEQUENCES FOR SGRBS AND 
X-RAY AFTERGLOWS 

The timescale over which the baryon-loaded wind and its 
associated EM luminosity can be sustained depends on the 
timescale for removal of differential rotation and is uncertain. 
An order of magnitude estimate is provided by the Alfven 
timescale m, which corresponds to ~ 1 — 10 s for the mag¬ 
netic fields considered here and which is in agreement with 


an estimate fl/Cl < 10 s for the timescale of change of the 
angular velocity profile in our simulations. Differential rota¬ 
tion is completely removed within a few to several of these 
timescales. The overall duration of the process as well as 
the luminosities mentioned above are compatible with the ob¬ 
served extended emission of SGRBs (lasting < 100 s) and 
thus make the magnetically-driven wind a promising alterna¬ 
tive model to explain such emission. 

The prompt SGRB emission is thought to be collimated and 
thus beamed away from us in the vast majority of events. Fur¬ 
thermore, the observed SGRB rate within the sensitivity vol¬ 
ume of advanced LIGO and Virgo is very small (< lyr -1 
ED- Therefore, identifying potential isotropic EM emission 
from BNS mergers is of utmost importance. The isotropic EM 
emission found here potentially represents such an important 
counterpart to the GW signal of BNS mergers. 

Although the initial phase of differential rotation of the 
BNS merger remnant is rather short, the total amount of mass 
lost via the magnetically-driven wind can be substantial. The 
ejected matter creates an optically thick environment that will 
likely affect the emission properties of the system in the fol¬ 
lowing evolution, once the star has settled down to uniform 
rotation. At this stage, the NS might lose rotational energy via 
spin-down radiation, and the resulting X-ray afterglow radia¬ 
tion would be the product of a reprocessing of this radiation 
by the confining envelope of ejecta matter (cf. flU ). The 
associated spectrum could carry such a signature, e.g., ther- 
malization of (part of) the non-thermal spin-down emission 
by the optically thick ejecta material. 

As the presence of a magnetically-driven wind only de¬ 
pends on differential rotation, which is generic to BNS merger 
remnants, the associated baryon pollution in the polar regions 
is also generally expected. This severely hampers the produc¬ 
tion of a relativistic jet through accretion of a remnant torus 
around the newly formed NS. Hence, if the presence of long- 
lasting X-ray afterglows requires a long-lived NS remnant, an 
alternative explanation for the generation of the SGRB itself 
is likely required. Such a scenario that explains the production 
of long-lasting X-ray afterglows and the prompt SGRB emis- 
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sion in a common phenomenology has recently been proposed 
fm (for an alternative proposal, see HE)). However, full sim¬ 
ulations of BNS mergers will be needed in the future to better 
characterize the magnetically-driven wind. These simulations 


will also help to shed light on whether or not a relativistic jet 
can be launched at the merger of a BNS system forming a 
long-lived NS. 
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